Enhancing piezoelectricity through polarization-strain coupling in ferroelectric 

superlattices 
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Short period ferroelectric/ferroelectric BaTiOs (BTO)/PbTi03 (PTO) superlattices are studied 
using density functional theory. Contrary to the trends in paraelectric/ferroelectric superlattices 
the polarization remains nearly constant for PTO concentrations below 50%. In addition, a signif- 
icant decrease in the c/a ratio below the PTO values were observed. Using a superlattice effective 
Hamiltonian we predict an enhancement in the 1^33 piezoelectric coefficient peaking at ~75% PTO 
concentration due to the different polarization-strain coupling in PTO and BTO layers. Further 
analysis reveals that these trends are bulk properties which are a consequence of the reduced P 
brought about by the polarization saturation in the BTO layers. 
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Perovskite superlattices present a new paradigm for en- 
gineering ferroelectrics and piezoelectrics for modern de- 
vice applications. Epitaxial strains, resulting from lattice 
mismatches, and induced changes in polarization, result- 
ing from electrostatic considerations, can have significant 
effects on the macroscopic properties of these artificial 
structures. 0, Q Superlattices combining a paraelectric 
(PE), such as SrTi0 3 (STO), with a ferroelectric (FE), 
such as PbTi0 3 (PTO) or BaTi0 3 (BTO) have been well 
studied both by experiment and theory. 0, 0, H, H, 
Their main features can be modeled by considering them 
as layers of strained bulk-like material with appropriate 
electrostatic boundary conditions in some systems, 
interface effects have also been shown to be important in 
the limit of ultrathin layers Q. In this paper, we extend 
this first-principles modeling approach to superlattices 
combining two ferroelectrics, which we expect to exhibit 
new features arising from the interplay of the different 
bulk spontaneous polarizations and dielectric and piezo- 
electric coefficients of the two constituents. 

We perform first-principles calculations, using density 
functional theory (DFT) for PTO/BTO superlat- 

tices with m=l-3 and n=l-3 layers of PTO and BTO, 
respectively. All calculations used projector augmented 
wave (PAW) potentials 0, lll| with the Vienna Ab ini- 
tio Simulation Package (VASP v4.6.26)) [H, witn the 
local density approximation for the exchange correlation 
functional. A 700 eV (22 Ha) cutoff and a 6x6xZ k- 
point mesh were used (where l=i for periods 2 and 3 
superlattices, 2 for periods 4 and 5 and 1 for period 
6). For BTO in the tetragonal P4mm 5-atom unit cell 
structure, the computed lattice constants are a=3. 947 A 
and c=3.996 A. The tetragonal PTO lattice constants 
were computed as a=3. 867 A and c=4.033 A [experi- 
ment: a=3.904 A and c=4.152 A].fl3| This agreement is 
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FIG. 1: Polarization along the z axis (top) and c/a ratio 
(bottom) as a function of % PTO for BTO/PTO superlat- 
tices with varying superlattice periods. Solid squares repre- 
sent DFT calculations and the dashed lines are the values 
obtained using the superlattice model (see Eqs. [1] and [3j . 
Numerical model parameters can be found in Table H] 



typical of LDA calculations for ferroelectric perovskites. 
In all superlattice calculations the in-plane lattice con- 
stant was constrained to that of the theoretical value 
for an STO substrate (3.863 A) while the c lattice vec- 
tors were optimized within the PAmm space group with 
lxl in-plane periodicity. Additional calculations for a 
doubled -^2x^2 in-plane unit cell for m=n=l show no 
instability to octahedral tilts of the type discussed in 
Ref. lZ|. All ionic coordinates were fully relaxed until 
the Hellman-Feynman forces on the ions were less than 
5 meV/A. Polarizations were computed using the Berry 
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phase method [141] . For bulk PTO, constrained to the 
hypothetical STO substrate in-plane lattice constant of 
3.863 A, the out-of-plane lattice constant was 4.039 A 
and the polarization was 0.79 C/m 2 . The epitaxial con- 
straints placed on the BTO layers result in a BTO c/a 
ratio (1.063) much larger than that of the more polar 
PTO layers. The BTO polarization was 0.39 C/m 2 , 
substantially enhanced over the computed bulk value of 
0.26 C/m 2 . 

Figure [T] shows the dependence of computed polariza- 
tion and c/a lattice constant on the concentration of PTO 
and the superlattice period of short period PTO/BTO 
superlattices for the full set of structures with period < 
6. Consistent with the fact that PTO has a greater po- 
larization than BTO, higher concentrations of PTO give 
rise to larger P of the superlattice. Similarly, decreasing 
the content of the larger c/a ratio strained BTO reduces 
the average c/a ratio of the superlattice. However, closer 
examination shows some unexpected features. First, the 
polarization of the superlattices remains almost constant 
for PTO concentrations below 50 % while the c/a ratio 
is below, rather than above, that of bulk PTO for a wide 
range of PTO concentrations (~40 % to 100%). In fact, 
the derivative of c/a with respect to PTO concentration 
is positive near 100%, rather than negative as would be 
expected from a simple linear interpolation. 

To interpret our DFT data and to distinguish between 
interface effects and bulk properties, we constructed a 
parameterized energy expression for the superlattice fol- 
lowing Refs. [HI and O. The expression is a Taylor 
expansion around the cubic perovskite structure in terms 
of the six independent components rji of the strain ten- 
sor (i is a Voigt index, i =1-6) and the three Cartesian 
soft mode amplitude components u a (a — x,y, z) param- 
eterized from DFT calculations. Restricting to PAmm 
symmetry, the effective stress-strain elastic enthalpy for 
a pure component ferroelectric under epitaxial strain (i.e. 
T}x = T)2 = ff) with out-of-plane stress, 173, can be written 
as: 



G(fj, 7)3, u z , a 3 ) = -Bii(2fy 2 + % 2 ) + B 12 (if + 2 m fj) + 



ku, + au, 



~Bi xx r) 3 u z + B lyy rju 2 z - (T3773 ,(1) 



where B\\ and B\ 2 are related to the elastic constants 
of the crystal, k and a are two independent symmetry- 
allowed fourth-order coefficients describing the cubic 
anisotropy and B\ xx and B\ yv are the phonon-strain cou- 
pling coefficients. 

Figure [2] shows the polarization versus in-plane lat- 
tice constant for pure component BTO and PTO. The 
model parameters (TableHJ were fit to DFT data for PTO 
with in-plane lattice constants within ±1% of the STO 
in-plane lattice constant 3.863A . For BTO, the fitting 
range was extended to include the computed bulk tetrag- 



O 
o 

"4— ' 

03 
N 

jjj 

O 
Q_ 



1.00r— 
0.90 f*- 
0.80 - 
0.70 - 
0.60 - 
0.50, r 
0.40 - ' 
0.30- 
0.20- 
0.10- 
0.00 1— 



DFT (PTO) 
DFT (BTO) 
Effective Model 



3.84 3.88 3.92 = 3.96 

In-plane constant [A] 



4.00 



FIG. 2: Polarization along the z axis as a function of in- 
plane strain for bulk BTO and PTO. Solid circles and squares 
represent DFT calculations for PTO and BTO, respectively. 
Dashed lines are the corresponding values obtained using Eq. 
[JJ with the parameters of Table |TJ 



onal lattice constant of 3.947A . The excellent agreement 
between the model and our DFT calculations is an indi- 
cation of the quality of the model. 

For a given mpTO /^bto superlattice we assume uniform 
soft mode amplitude u z within each constituent layer. 
Further, we impose the constraint that the correspond- 
ing polarizations of the two constituent layers i = BTO, 
PTO, 



PI = —Z*' l ul 
z W 



(2) 



are equal, where e is the absolute value of the electron 
charge, fi l is the layer unit cell volume and Z*' 1 is the 
Born effective charge of the soft mode, so that polar- 
ization is uniform throughout the superlattice. Strictly 
speaking, it is the displacement field, D = P + sE, which 
remains uniform throughout the superlattice. How- 
ever, this uniform polarization approximation has pre- 
viously been shown to be valid for short-period super- 
lattices [H, [TtJ and is supported in the present case by 
unit-ccll-laycr polarization profiles computed using bulk 
Born effective charges (not shown). 

The effective elastic enthalpy for the two-component- 
superlattice with mpTO layers and tibto layers is then 
obtained as: 



TABLE I: Energy expansion coefficients, Eq.fJJ and soft mode 
Born effective charges, for PTO and BTO, in atomic units fit 
to the DFT results in Fig. [2] as described in the text. 
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do 


BaTiOs 


5.13 


3.05 


-1.00 


-0.100 


-0.007 


0.14 9.94 


7.48 


PbTiOs 


4.69 


1.15 


-0.705 


0.207 


-0.0132 


0.0364 9.40 


7.385 
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FIG. 3: Dependence of a?33 on PTO concentration for 
PTO/BTO superlattices. The dashed line represent the su- 
perlattice effective Hamiltonian at bulk STO in plane lattice 
constant of 3.863 A and the dotted line is are the same for 
an in-plane lattice constant of 3.90 A. The solid squares are 
computed from DFT. 
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(3) 



The superlattice P and 7? 3 are determined by minimiz- 
ing Eq.[3]with respect to 773 and u l z , imposing the uniform 
polarization approximation through the relation 



,BTO 



rwP TO ,withr 



n BTO z 



*.PTO 



^PTO^* : BTO ' 



(4) 



This model can also be used to compute the piezoelec- 
tric coefficient, CZ33, as: 



dP 
da 3 



meZ*< PTO B[Z° 



2n PTO u PTO (2B[ 1 TO a PTO 



meZ 



*,BTO 



BTO 

xx 



2n BTO uf TO (2B PTO a BTO - \{B{ 



}BTO\ 
lxx 1 



(5) 



Figure [3] shows the dependence of g?33 on composition as 
predicted from Eq. O The model is in good agreement 
with the DFT d 33 values for bulk BTO (32 pC/N), bulk 
PTO (55 pC/N) and the 3 PTO / 1 BTO superlattice 
(58 pC/N). Surprisingly, we find a 5% enhancement of the 
e?33 coefficient at 75 % PTO where the c/a is a minimum. 

The good agreement between the DFT results and the 
superlattice energy functional, shown in Fig. [TJ suggests 
that the observed trends in P and c/a are a consequence 
of bulk electrostatics and strain effects, and that inter- 
faces do not play an essential role. The success of the 
model allows us to elucidate the origin of the nearly con- 
stant P for low PTO concentration through an exami- 
nation of the evolution of the FE potential energy well 



(Fig. [H left panel) . A comparison of the FE potential en- 
ergy wells of the epitaxially strained parent compounds 
shows that BTO has a much stiffer FE well with a mini- 
mum at a P of less than half of PTO. For polarizations 
much greater than 0.39 C/m 2 the effective enthalpy of 
BTO sharply increases, while PTO has a much gentler 
dependence. Here the P in the BTO layers is saturated, 
requiring larger electric fields to further polarize them. 
Similar polarization saturation effects have been observed 
in compressively strained PZT and PTO L 18j|. In a su- 
perlattice, this translates into the need for a higher per- 
centage of the more polar PTO in order to increase the 
macroscopic polarization of the superlattice. This fact is 
borne out in the evolution of the FE wells with % PTO. 
Even up to PTO concentrations of 80% the FE well still 
resembles that of BTO, severely limiting the total P in 
the superlattice and accounting for the suppression in 
expected P enhancement with increasing % PTO. 

This suppressed layer-by-layer polarization can be 
linked to the abnormal decrease of c/a below that of 
the PTO parent compound (Fig. |4|) as well as the en- 
hancement in the d 33 coefficient. The model indicates 
that if P were to increase linearly with increasing PTO 
concentration then the c lattice parameter would have a 
corresponding linear decrease to the PTO values. How- 
ever, as previously stated, the P throughout the super- 
lattices is severely suppressed due to the presence of the 
BTO layers. Since P is coupled to the c lattice parameter 
this results in a drastic reduction in the average c lattice 
parameter (see Fig. [4]). This BTO effect is only over- 
come at extremely high concentrations of PTO, where 
the electric field generated by the PTO layers are capa- 
ble of counteracting the BTO potential energy well con- 
straints. Furthermore, BTO has a much steeper depen- 
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FIG. 4: Effective enthalpy (left) and c lattice parameter 
(right) as a function of polarization for various compositions 
of PTO /BTO superlattices obtained from the superlattice en- 
ergy functional. Dashed lines represent pure BTO, dashed- 
dotted lines are for pure PTO and solid lines indicate the 
ferroelectric wells for the superlattices at 20% PTO intervals. 
Open circles mark the c lattice parameter for each composi- 
tion at the predicted P. The dotted line in the right figure 
indicates c for epitaxially-strained tetragonal PTO. Energies 
are relative to epitaxially strained bulk BTO. 
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dence of c on P than PTO, i.e. the c lattice parameter 
increases more with increasing P in the BTO layers than 
in the PTO layers. As such, it is the balance between the 
stronger polarization-strain coupling in the BTO layers 
and the amount of PTO present to pole these layers which 
give rise to the peak in ^33 at ~75% PTO concentration. 
These results suggest that it may be possible to further 
enhance by reducing the saturation in the BTO layers 
either by decreasing the in-plane compressive strain (see 
Fig. IU dotted line: a=3.90A) or through other factors 
such as intermixing at the interface between PTO and 
BTO layers. @ 

In conclusion, we have used density functional theory 
to explore the dependence of P and c/a as a function 
of PTO concentration for short period PTO/BTO su- 
perlattices. Our extensive DFT calculations reveal two 
intriguing phenomena: nearly constant P for PTO con- 
centrations less than 50% and a dramatic decrease in c/a 
to values below that of the smaller PTO. Using a super- 
lattice effective Hamiltonian we demonstrate that these 
trends are a consequence of P saturation in the BTO 
layers which, by limiting the magnitude of P, results in 
abnormally large decreases in the c lattice parameter. 
The competition between the polarization saturation ef- 
fects of BTO and the high polarization in the PTO lay- 
ers results in a corresponding peak in the coefficients 
in the PTO/BTO superlattices, suggesting a new way 
of enhancing the piezoelectric properties of the super- 
lattices; other material combinations may yet be found 
which show an even greater enhancement based on this 
mechanism. In addition, this model can be easily ex- 
tended to multicomponent systems, different phases and 
to model strain effects on these phases; allowing for the 
easy exploration of materials properties. 
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